Bifurcation in ground-state fidelity and universal order parameter for two-dimensional quantum 

transverse Ising model 
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We establish an intriguing connection between quantum phase transitions and bifurcations in the ground-state 
fidelity per lattice site, and construct the universal order parameter for quantum Ising model in a transverse 
magnetic field on an infinite-size square lattice in two spatial dimensions, a prototypical model with symmetry 
breaking order. This is achieved by computing ground-state wave functions in the context of the tensor network 
algorithm based on the infinite projected entangled-pair state representation. Our finding is applicable to any 
systems with symmetry breaking order, as a result of the fact that, in the conventional Landau-Ginzburg- Wilson 
paradigm, a quantum system undergoing a phase transition is characterized in terms of spontaneous symmetry 
breaking captured by a local order parameter. In addition, a bifurcation in the reduced fidelity between two 
different reduced density matrices is also discussed. 
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I. INTRODUCTION 

In recent years, a novel approach to quantum phase tran- 
sitions (PQTs) HI 01 in quantum many-body lattice systems 
emerges, based on the fact that fidelity, a basic notion in quan- 
tum information science, is a measure of quantum state distin- 
guishability As argued in Refs. , the ground-state 

fidelity per lattice site captures drastic changes of the ground- 
state wave functions around a critical point for a quantum sys- 
tem in any spatial dimensions. That is, it is a universal marker 
to detect QPTs, regardless of what type of the internal order 
is present. In fact, quantum many-body lattice systems with 
either spontaneous symmetry breaking (SSB) order or topo- 
logical order, have been specifically studied: for the former, 
the examples include quantum transverse Ising ||8|] and Potts 
chains IU3I1 . quantum Ising ]7[] and XYX models IU4I1 on a 
square lattice, while for the latter, the Kitaev model IU5I1 on a 
honeycomb lattice. Remarkably, there is a smoking-gun sig- 
nature for SSB order in the fidelity per site approach: SSB 
order implies a bifurcation, arising from degenerate ground- 
state wave functions due to broken symmetry, in the ground- 
state fidelity per lattice site. 

In fact, it was shown that, for a QPT arising from an SSB, 
a bifurcation appears in the ground-state fidelity per site, with 
a critical point identified as a bifurcation point. This in turn 
results in a novel concept of the universal order parameter, 
which appears as the ground-state fidelity per lattice site be- 
tween a ground state and its symmetry-transformed counter- 
part. The advantage of the ground-state fidelity per lattice site 
and the universal order parameter [16] over local order pa- 
rameters lies in the fact that both of them are universal, in 
the sense that it is not model-dependent, in contrast to model- 
dependant order parameters in characterizing QPTs in quan- 
tum lattice many-body systems. Similarly, a bifurcation oc- 
curs in the ground-state reduced fidelity [17] between the one- 
site reduced density matrices and the two-site reduced den- 
sity matrices, with a critical point identified as a bifurcation 
point. However, all specific examples, up to now, have been 



restricted to quantum many -body lattice systems in one spatial 
dimension. 

In this paper, we take one step further to see if these novel 
ideas are practical to quantum many-body lattice systems in 
two spatial dimensions. We establish, by exploiting the tensor 
network algorithm based on the infinite projected entangled- 
pair states (iPEPS) algorithm 111 811 . an intriguing connection 
between a QPT and a bifurcation in the ground-state fidelity 
per lattice site, and construct the universal order parameter for 
quantum transverse Ising model on a square lattice, a proto- 
typical model with symmetry breaking order. Our finding is 
applicable to any systems with symmetry breaking order, as a 
result of the fact that, in the conventional Landau-Ginzburg- 
Wilson paradigm, a quantum system undergoing a phase tran- 
sition is characterized in terms of SSB captured by a local 
order parameter. In addition, a bifurcation in the reduced fi - 
delity between two different reduced density matrices is also 
discussed. 



II. MODEL 

We consider quantum Ising model in a transverse magnetic 
field on an infinite-size square lattice in two spatial dimen- 
sions. It is described by the Hamiltonian: 



(i) 



(if) 



where cr' a (a = x,z) are the spin 1/2 Pauli operators at site 
i, (ij) runs over all the possible nearest-neighbor pairs on a 
square lattice, and A is the transverse magnetic field, which we 
choose as a control parameter. The model is invariant under 
the symmetry operation: ay — » -ay and cr'' 1 — > o4 !l for 
all sites simultaneously, which yields the Z2 symmetry. As is 
well known, the system undergoes a second-order QPTs at the 
critical field A c ~ 3.044 111 El. 
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III. BIFURCATION IN THE GROUND-STATE FIDELITY 
PER LATTICE SITE 

For two ground-state wave functions ^(Ai)) and 
|*F(/i2)) corresponding to two different values A\ and 
A 2 of the control parameter A, the ground-state fidelity 
F(A U A 2 ) = K < PW2)l 4 I'(^i)>l asymptotically scales as 
F{A\,A2) ~ d{A\,A2) N , with N being the number of sites 
in the lattice. Here, d{A\,A 2 ) is the scaling parameter, 
introduced for one-dimensional quantum systems |4j-|6J] and 
for two and higher-dimensional quantum systems H. Note 
that it characterizes how fast fidelity goes to zero when the 
thermodynamic limit is approached. Physically, the scaling 
parameter d(A\,A?) is the averaged fidelity per lattice site, 
which is well defined in the thermodynamic limit: 



i.oFTTTrrri: — 
*•••• 



In d(A\,A 2 ) = lim 

N— >oc 



In F(A U A 2 ) 
N ' 



(2) 



As noted in Refs. |4j-|7[], it satisfies the properties inherited 
from the fidelity F(A\,A 2 ): (i) normalization d(A,A) = 1; 
(ii) symmetry d(Ai,A 2 ) = d(A 2 ,A\)\ and (iii) range < 
d{A x ,A 2 )<\. 

In the Z 2 symmetric phase, the ground-state wave function 
is non-degenerate, while in the Z 2 symmetry-broken phase, 
two degenerate ground-state wave functions arise. If we 
choose ^¥{A 2 )) as a reference state, with A 2 in the Z 2 sym- 
metric phase, then the ground-state fidelity per lattice site, 
d(Ai,A 2 ), is not able to distinguish two degenerate ground- 
state wave functions in the Z 2 symmetry-broken phase. How- 
ever, if we choose PF^)) as a reference state, with A 2 in the 
Z 2 symmetry-broken phase, then the ground-state fidelity per 
lattice site, d(Ai , A 2 ), can be used to distinguish two degen- 
erate ground-state wave functions. Therefore, a bifurcation 
occurs in the ground-state fidelity per lattice site, d(A\,A2), as 
a function of A\, for a fixed A 2 . As argued in Refs. [8, 9], 
for a given truncation dimension D, a pseudo-phase-transition 
point A-q manifests itself as a bifurcation point 112011 . 

In Fig. [U we plot the ground-state fidelity per lattice site, 
d(Ai,A 2 ), for quantum Ising model in a transverse magnetic 
field on an infinite-size square lattice in two spatial dimen- 
sions, with the transverse magnetic field A as the control pa- 
rameter. Note that a bifurcation does occur for the ground- 
state fidelity per lattice site, d{A\,A 2 ). Here, we have chosen 
Ai = 2.1 as a specific example. The pseudo-phase-transition 
point locates at Ad - 3.100 for the truncation dimension 
D = 2, at Aj} = 3.065 for the truncation dimension D = 3, 
and at Ajj> - 3.050 for the truncation dimension D = 4. We 
mention that the ground-state fidelity per lattice site, d(A\ , A 2 ), 
is computed from the iPEPS representation of the ground- 
state wave functions, following the transfer matrix approach 
described in Ref.j7|]. 

A remarkable feature of the bifurcation points for the 
ground-state fidelity per lattice site, d{A\ , A 2 ), as seen in Fig.Q] 
is that d{A\,A 2 ) between different symmetry breaking ground- 
state wave functions in the same phase appears less than 
that between two ground-state wave functions from different 
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FIG. 1: (color online) The ground-state fidelity per lattice site, 
d(Ai,A.2), for quantum Ising model in a transverse field on a square 
lattice in two spatial dimensions, with the transverse magnetic field 
A as the control parameter. If we choose | v P(/t 2 )> as a reference state, 
with A 2 being in the Z 2 symmetry-broken phase, then d(A[ , A 2 ) distin- 
guishes two degenerate ground-state wave functions, with a pseudo- 
phase-transition point A B as a bifurcation point. Here, we have cho- 
sen A2 = 2.1. The pseudo-phase-transition point is identified at 
/Id = 3.100 for the truncation dimension D = 2, at A B = 3.065 for 
the truncation dimension D = 3, and at Ad = 3.050 for the truncation 
dimension D = 4, respectively. 



phases. This is due to the fact that two degenerate symmetry 
breaking ground states in the same symmetry-broken phase 
are more distinguishable than two ground states from differ- 
ent phases. 



IV. THE GROUND-STATE REDUCED FIDELITY BETWEEN 
TWO REDUCED DENSITY MATRICES 

For quantum Ising model in a transverse magnetic field on 
an infinite-size square lattice, the one-site reduced density ma- 
trix in the Z 2 symmetric phase takes the form, 



p=^+2(S z )S z , 



(3) 



where (S z ) is the ground-state expectation value of S z in the 
Z 2 symmetric phase, while the two-site reduced density matrix 
is, 



+ 4y xx S x ®S X + Ay zz S z ® S z 
+y oz I®S z +y zo S z ®I. 



(4) 



Here, y xx = {S x ® S x ), y zz =(S Z ® S z ), y oz = </ ® S z ), y zo = 
(S z <8> /), and / is the identity matrix. 

In the Z 2 symmetry -broken phase, the one-site reduced den- 
sity matrix becomes, 



p = ^+2(S x )S x + 2(S t )S z , 
whereas the two-site reduced density matrix is 

p = ^7 + 4y xx S x ®S x +4y zz S z ®S z 

+y oz I ®S Z + y zo S Z ®I + 4y xz S X ®S Z 
+4y zx S Z ®S X + y ox l ®S X + y xo S X ®I, 



(5) 



(6) 
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with y xz = (S x ® S z >, y ZA = (S : <g> S^), y ox = <7 ® S A ->, and 

r.to = (s x ®r>- 

The reduced fidelity measures the distance between two 
quantum mixed states. For two reduced density matrices p A , 
and p A „ the reduced fidelity F(p Ai ,p A2 ) is defined to be 11711 

F(pa i ,Pa 2 ) = tryjp^p^p^ 2 . (7) 

Here, p^, and p Al are the reduced density matrices correspond- 
ing to two different values, A\ and Xi, of the control parameter 
A. Notice that the reduced fidelity F{p^,p^) is a function of 
/I] and /I2, which satisfies the following properties: (i) normal- 
ization F(p /U p A ) = 1; (ii) symmetry F(p M ,p Al ) = F(p Al ,p Al ); 
(hi) range 0< F(p, h ,p Al ) <1. 

In Fig. ?? (upper panel), we plot the ground-state reduced 
fidelity F(p Al ,p Al ) between the one-site reduced density ma- 
trices for quantum Ising model in a transverse field on an 
infinite-size square lattice, with the transverse field strength 
A as the control parameter. Here, we choose p Al , at Ai = 2.1, 
as a reference state, which breaks the Z2 symmetry. The one- 
site reduced fidelity is able to distinguish two mixed states 
(described by two reduced density matrices) from two degen- 
erate ground-state wave functions, with a bifurcation point 
as a pseudo-phase-transition point Ad- When the control pa- 
rameter A crosses a pseudo-transition point, the ground-state 
degeneracy changes suddenly, implying that the system un- 
dergoes a QPT We observe that the pseudo-phase-transition 
point Ac moves toward the critical point A c with increasing D. 
More precisely, the pseudo-phase-transition point locates at 
Ad = 3.100 for the truncation dimension D = 2, at Ajy = 3.065 
for the truncation dimension D = 3, and at Ad = 3.050 for the 
truncation dimension D = 4, respectively. 

We also plot the ground-state reduced fidelity F(p Al ,p Al ) 
between the two-site reduced density matrices for quantum 
Ising model in a transverse field on an infinite-size square lat- 
tice in Fig. ?? (lower panel). The reference state is chosen 
at the same A2 = 2.1, as in the case of the one-site reduced 
fidelity. We observe that a bifurcation also occurs in the two- 
site reduced fidelity, with the same pseudo-phase-transition 
point A®. This is expected, simply because they are resulted 
from the same set of the ground-state wave functions. 

V. THE UNIVERSAL ORDER PARAMETER 

As argued in Ref. lfl6ll . for any quantum lattice system with 
a symmetry group G undergoing a QPT with symmetry break- 
ing order, there is a universal order parameter: it is defined 
as the ground state fidelity per lattice site between a ground- 
state wave function and its symmetry-transformed counter- 
part, which is discontinuous for first-order phase transitions 
and continuous for second-order phase transitions. This is 
based on the observation that, for any ground state IT*) in the 
symmetric phase, ( v F|g| v F) is equal to 1, for any symmetry op- 
eration g e G, whereas it is identical to zero for any state in 
the symmetry broken phase. 




FIG. 2: (color online) Upper panel: The ground-state one-site re- 
duced fidelity, F(p Ai ,p Al ), for quantum Ising model in a transverse 
field on a square lattice in two spatial dimensions. Here, we have 
chosen p Al as the reference state, with A2 = 2. 1 in the Zi symmetry- 
broken phase. Lower panel: The ground-state two-site reduced fi- 
delity, F(p Al ,p Al ), for quantum Ising model in a transverse field on 
a square lattice. The reference state p Al has been chosen at the same 
A2 = 2.1 as for the one-site reduced fidelity. Note that F(p Al ,Pa 2 ) is 
able to distinguish two mixed states arising from degenerate ground- 
state wave functions. A pseudo-phase-transition point occurs as a 
bifurcation point. In both cases, the pseudo-phase-transition point 
locates at /t D = 3.100 for the truncation dimension D = 2, at 
Ad = 3.065 for the truncation dimension D = 3, and at Ad = 3.050 
for the truncation dimension D = 4, respectively. 

In order to measure the distance between two quantum 
states and gf¥(X)), let us consider their counterparts 

|T(A))/, and g^iA))^ on a finite-size lattice, with L being 
the number of the total lattice sites. As argued in Refs. [4- 
7], Il( v P|§I*1')lI asymptotically scales as fHX) with L, as one 
may see from the tensor network representations of the sys- 
tem's ground state wave functions. Here, f g (A) is the aver- 
aged fidelity per lattice site, which is well-defined even in 
the thermodynamic limit. As such, one sees that, f g (A) = 1 
for any g s G, if A is in the symmetric phase A > A c , and 
< f g (A) < l for any nontrivial symmetry operation g, if A is 
in the symmetry -broken phase A < A c . As argued in Ref. itlol . 
we define the universal order parameter to be 

I g (A) = ^Jl-fi(A). (8) 

Note that I g {A) is always zero if A > A c . However, it becomes 
nonzero, with its value ranging from to 1, if A < A c . These 
features are exactly what one requires for I g (A) to be an order 
parameter. In fact, this is valid for any quantum many-body 
lattice system with a global symmetry group G spontaneously 
broken. A remarkable feature of the universal order parameter 
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FIG. 3: (color online) The universal order parameter 1(A) for quan- 
tum Ising model in a transverse magnetic field on a square lattice. 
A pseudo-phase-transition point Ac occurs, as 1(A) changes from be- 
ing nonzero to zero at A = A®. When the truncation dimension D is 
increased, a pseudo-phase-transition point A B approaches the critical 
point A c . In addition, the universal order parameter 1(A) reaches the 
maximum value 1(A) = 1 at the factorizing field Af = 0. 



for quantum Ising model in a transverse magnetic field on an 
infinite-size square lattice. We expect that our approach might 
provide further insights into critical phenomena in quantum 
many-body lattice systems in condensed matter. 
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is that it not only makes it possible to locate a critical point, 
but also enables us to identify a factorized state ^(Af)), with 
Af being the so-called factorizing field 112 ill . 

In Fig[3] we plot the universal order parameter I g (A) for 
quantum Ising model in a transverse field on an infinite-size 
square lattice, with the field strength A as the control param- 
eter. Here, the symmetry operation is the non-trivial element 
g of the group Z2. If A < Ad, the universal order parameter 
I g (A) is non-zero. This characterizes the Z2 symmetry-broken 
phase, in contrast to the fact that the universal order parame- 
ter I g (A) is zero, in the symmetric phase A > Ad- When the 
control parameter A varies across the pseudo-critical point Aq, 
the behavior of the universal order parameter I g (A) changes 
qualitatively, implying that the system undergoes a QPT at 
the pseudo-phase-transition point Aq. As the truncation di- 
mension D is increased, the pseudo-phase-transition point Ajy 
moves toward the critical point A c . In addition, the (trivial) 
factorizing field Af = exists for quantum Ising model in 
a transverse field on an infinite-size square lattice, at which 
I g (A) reaches its maximum. 

VI. CONCLUSIONS 

We have investigated an intriguing connection between 
QPTs and bifurcations in the ground-state fidelity per lattice 
site, in the context of the tensor network algorithm based on 
the iPEPS representation. For quantum transverse Ising model 
on an infinite-size square lattice, the iPEPS algorithm pro- 
duces two degenerate symmetry-breaking ground-state wave 
functions arising from the Z2 symmetry breaking, each of 
which results from a randomly chosen initial state. There- 
fore, a quantum system undergoing a phase transition is char- 
acterized in terms of SSB that is captured by a bifurcation in 
the ground-state fidelity per lattice site. We have also con- 
structed the universal order parameter , and discussed bifurca- 
tions in the ground-state reduced fidelity between two differ- 
ent reduced density matrices in the symmetry-broken phase, 
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